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ABSTRACT 

Using the Planck far-infrared and Arecibo GALFA 21-crn line surveys, we identified 
a set of isolated interstellar clouds (approximately degree-sized on the sky and com¬ 
prising 100 solar masses) and assessed the ratio of gas mass to dust mass. Significant 
variations of the gas-to-dust ratio are found both from cloud to cloud and within regions 
of individual clouds; within the clouds, the atomic gas per unit dust decreases by more 
than a factor of 3 compared to the standard gas-to-dust ratio. Three hypotheses are 
considered. First, the apparently low gas-to-dust ratio could be due to molecular gas. 
Comparing to Planck CO maps, the brightest clouds have a H 2 /CO ratio comparable 
to galactic plane clouds, but a strong lower limit is placed on the ratio for other clouds, 
such that the required amount of molecular gas is far higher than would be expected 
based on the CO upper limits. Second, we consider self-absorbed 21-crn lines and find 
the optical depth must be ~ 3, significantly higher than found from surveys of radio 
sources. Third, grain properties may change within the clouds: they become more emis¬ 
sive when they are colder, while not utilizing heavy elements that already have their 
cosmic abundance fully locked into grains. It is possible all three processes are active, 
and follow-up studies will be required to disentangle them and measure the true total 
gas and dust content of interstellar clouds. 

Subject headings: dust, ISM: abundances, ISM: atoms, ISM: clouds, ISM: general, 

ISM: molecules 



1. Introduction 


The nature of the diffuse interstellar medium, which pervades the volume of the galaxy and is 
distinct from the dense, self-gravitating clouds where stars form, is only understood in broad terms. 
Studies of absorption toward nearby stars reveal a combination of atomic and molecular gas toward 
those lines of sight ( Savage et al. 1 1977 Rachford et al.|2002 ), but they leave large gaps between the 
observed stars and have relatively little context by which we can understand the organization of the 
material or its relation to material not readily studied by absorption lines. Large-scale surveys of 
emission from diffuse interstellar gas were made possible by the brightness of the 21-cm hyperfine 


line of H I and dedicated all-sky surveys at half-degree angular resolution (Kalberla et al. 2005) 


Widespread emission from interstellar dust was revealed by the first infrared satellite observations 
and unveiled even more detail about the organization of the interstellar medium (Low et al.|l984). 
The structure of the interstellar clouds and the relation between gas and dust could be studied 
at half-degree resolution, limited by the 21-cm surveys. Dust and gas well correlated at that 
scale; however, the infrared images showed significant structure on smaller scales, and variations 
between dust and gas appeared significant on small scales in limited observations performed. As one 
observer noted, “[o]ur results emphasize that it is not always evident, from observations with spatial 
resolutions coarser than a few arc minutes, how interstellar gas and dust are related” (Verschuur 


et al. 1992). 


At the level of individual clouds, where a ‘cloud’ is defined as a structure on the sky of high 
spatial contrast on the < 1° scale as well as a velocity coherent at the < 10 km s -1 scale, there are 
significant variations between the dust and gas properties. The advantage of working on features 
of relatively high contrast is that we can isolate the structures by subtracting the brightness just 
outside the cloud, allowing measurements of the cloud without requiring knowledge of the absolute 
goal was brightness. For degree-sized clouds selected from the Infrared Astronomical Satellite 
(IRAS), we found significant variations in the relative amounts of dust and gas from cloud to 


cloud (Heiles et al. 1988 

) despite the overall linear correlation between dust and gas when averaged 

over very large areas ( 

3oulanger & Perault 

1988 

). In particular, there was more dust emission 


per unit H I column density for some clouds, suggesting the dust was tracing additional material 
not traced by the 21-cm line observations. Follow-up observations showed CO in several of those 
clouds indicated the extra infrared emission was often associated with molecular gas, though the 
correspondence was not perfect ( Reach et al.|1994 Meyerdierks & Heithausen 1996). Some galactic 
CO emission is apparently present in small clumps, even at high galactic latitude (Heithausen 2002). 


All-sky comparisons of dust and gas emission, include far-infrared observations out to 240 
fim wavelength with the Cosmic Background Explorer at degree scales showed that locations with 
apparent ‘excess’ infrared emission tended to also be colder than their surroundings (so the excess is 
truly in the amount of dust, not the amount of heating) and provided an all-sky guide to locations 
of potential molecular clouds (Reach et al.| |l998). More recently, the all-sky survey by Planck with 
its exceptionally high sensitivity at 5' angular scales, confirmed the presence and relative frigidity 
of the infrared excess and showed that the amount of material that is implicated in producing 
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the excess is comparable to the known molecular mass of the local interstellar medium (Planck 
Collaboration XIX 2011). 


The potential for molecular gas without CO has been noted for some time (Blitz et al. 1990), 
but the true amount and extent of such gas is only recently becoming elucidated by observations 
such as mentioned above and described in this paper. H 2 is self-shielded and forms in regions 
protected front he interstellar radiation field by approximately half a magnitude of extinction, 
while CO requires somewhat more protection at approximately 2 magnitudes of extinction. A 
theoretical prediction for the amount of CO-dark molecular gas yields results that can generally 
agree with those inferred from observations, with a reasonable set of physical properties of the 
interstellar medium, depending primarily on the fraction of the total medium that is in very large 
clouds (Wolfire et al. 2010| ). Therefore, in regions with no giant molecular clouds, the fraction of 
‘dark’ molecular gas is expected to be high. The potential existence of untraced molecular gas 
is important for understanding the formation of molecular clouds on galactic scales, where large- 
scale molecular (but low-density) filaments are predicted to extend well beyond the observed dense 
filaments (Smith et al. 2014). 


In this paper, we study isolated structures in the local interstellar medium in order to compare 
dust and gas at 5' angular scales. We address the questions of whether dust and gas are well mixed, 
what is the relative amount of atomic and molecular gas in the interstellar medium, and where does 
the gas transform from atomic to molecular. 


Observations 


2.1. Arecibo 


The 21-cm line data used in this work were part of the large-scale Galactic Arecibo L-band 
Feed Array (GALFA) HI survey, Data Release 1 (Peek et al. 2011). The observations with the 305- 
m Arecibo telescope are part of a survey that will cover the entire sky visible from that telescope 
(declination 0° to 40°) with an angular resolution of 4'. The time-ordered spectra were regridded 
into cubes and presented in an online archive. We extracted cubes centered on distinct interstellar 
clouds evident in the Planck all-sky survey. An interference spike was removed from the cube near 
G94-36, and linear baselines were refit using symmetrically located portions on either side of the 
low-velocity H I that dominates the spectra. For each position, we integrate the 21-cm line and 
convert to a column density assuming that the line is optically thin. We return to address this 
assumption later in the paper. Note that the Arecibo observations contain stray-light due to 21-cm 
emission from distant sources (in particular, the galactic plane) entering the far side lobes of the 
beams of the 7 receivers. Comparison with the stray-light-corrected all-sky survey at low angular 
resolution (Kalberla et al. 2005) show the amplitude of stray light is less than ~ 1 K brightness 


temperature, which is much less than the brightness of the clouds studied in this paper (> 10 K). 
The stray-light should be smooth on the ~degree angular scale and smaller that are of interest in 
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this project. The stray-light should also have a relatively broad velocity profile (because it contains 
large portions of the rotating galactic plane) which we verify is not evident at the levels > 10 20 
cnr 2 of interest to this project. 


2.2. Planck 


The Planck survey covered the entire sky in 9 bands with wavelengths from 30 GHz to 857 GHz 
(1 cm to 350 /iin) ( Tauber et al.pOlO Planck Collaboration IpOll ). The wide range of frequencies 
allowed a separation of the sky into the cosmic microwave background and astrophysical foregrounds 
(Planck Collaboration XII 2014); we utilized the public release of the ‘thermal dust’ component. 
The optical depth of interstellar dust at 353 GHz was derived from a fit to the combination of 
the Planck 353, 545, and 857 GHz (wavelengths 850, 550, and 350 /im) from the High-Frequency 


Instrument (Planck Collaboration VIII 2014) together with the IRAS 100 fi m data as reprocessed 


to combine CORF/DIRBE large-scale calibration (Miville-Deschenes & Lagache 2005). After the 


contaminating signals from the CMB and zodiacal light, the dust emission was fitted with 


I u = t(353 GHz) 


353 GHz 


B v (T d ) 


( 1 ) 


where r is the optical depth and /3 is an emissivity power-law index that approximates the frequency- 
dependence of the dust emissivity. The fit for /3 was performed at 30' resolution, and then used to 
make a finer-scaled fit for T<i at the full 5' resolution. 


2.3. Cloud selection 


Our goal was to identify isolated, degree-sized clouds for which the unrelated background 
ISM could be separated from a distinct spatially compact (at 5' resolution) structure. The clouds 
studied in this paper were selected by inspecting the Planck 353 GHz map, in equatorial coordinates, 
within the Arecibo declination range — 1° < <5 < 38° and clearly separated from the ridge of emission 
centered in the galactic plane. Figure [I] shows the clouds on the Planck image, and Table[l]lists their 
locations. Despite being the prominent, isolated features in the Arecibo-visible portion of the sky, 
most of the clouds have not been previously studied. G94-36 includes MBM 53 and 54, which had 
been identified in a survey of the Palomar plates for visible extinction (Magnani et al. 1985). G104- 
39 includes DIR 105-38 , an infrared excess cloud identified from comparison of Cosmic Background 
Explorer /Diffuse Background Experiment (DIRBE) far-infrared observations to large-scale 21-crn 
survey (Reach et al.| 1998). 



















Table 1: Properties of the Interstellar Clouds in this Study 3, 


Name 

Coordinates 

(J2000) 

vo 

(km s” 1 ) 

SV 

(km s” 1 ) 

T d 

(K) 

P 

A1V(H I) 
(10 2 ° cm” 2 ) 

G94-36 

23:05:23 

+20:59 

-7.0 

4.90 + 0.08 

17.0 + 0.3 

1.78 + 0.01 

4.2 

G96-51 

23:43:14 

+08:58 

-9.1 

6.61 + 0.11 

18.0 + 0.3 

1.64 + 0.02 

4.3 

G104-39 

23:48:42 

+21:50 

-13.6 

6.8 ±0.2 

17.4 + 0.3 

1.79 + 0.02 

3.0 

G108-53 

00:17:06 

+09:45 

-5.7 

8.70 + 0.14 

17.0 + 0.2 

1.84 + 0.01 

4.6 

G138-52 

01:29:26 

+ 10:10 

-9.5 

6.45 + 0.05 

18.7 + 0.2 

1.71 + 0.01 

4.4 

G198+32 

08:27:26 

+26:12 

+5.1 

2.91 + 0.02 

18.4 + 0.7 

1.45 + 0.03 

3.0 

G198+32 E b 

08:16:15 

+27:02 

+5.1 

2.9 + 0.02 

19.3 + 0.3 

1.70 + 0.02 

1.5 

G221+35 

09:08:42 

+09:26 

+6.1 

10.1 + 0.2 

18.6 + 0.2 

1.66 + 0.01 

5.0 

G236+39 

09:45:58 

+00:35 

+10.2 

8.1 ±0.5 

18.0 + 0.3 

1.64 + 0.02 

3.0 

G236+39 N c 

09:46:46 

+01:51 

+9.3 

5.0 ±0.2 

18.9 + 0.2 

1.73 + 0.02 

5.7 

G254+63 

11:29:18 

+07:34 

-1.6 

6.22 + 0.09 

18.9 + 0.4 

1.25 + 0.01 

4.7 


a Td and /3 are the average dust temperature (K) and emissivity index over the cloud peaks, and the numbers after 
the ± are standard deviations within the cloud. AN(H I) is the 21-cm column density (in the optically-thin limit) 
toward the cloud peak minus a background within the 8° region containing the cloud. 

6 cloud near G198+32 with high gas/dust 
“northern component of G236+39 
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Fig. 1. - All-sky map of the Planck 353 GHz optical depth, r(353GHz), in celestial coordinates. 
The declination range not visible from Arecibo is greyed out. Each of the relatively isolated inter¬ 
stellar clouds addressed in this study are labeled with their name (based on galactic coordinates). 










Table 2: Derived properties 3 


Name 

?diff 

(K) 

\G/D] m 

-^cloud 

(K) 

[G/Oldoud 

/dark 

Cloud Mass 

(Mo) 

H I Mass 
(M Q ) 

G94-36 

18.5 ±0.6 

48 ± 11 

17.1 ±0.2 

17 ±2 

2.9 

106 

22 

G96-51 

18.3 ±0.4 

51 ±3 

17.5 ±0.2 

24 ± 3 

2.1 

72 

17 

G104-39 

18.4 ±0.5 

53 ±9 

17.4 ±0.2 

20 ±2 

2.6 

91 

14 

G108-53 

18.1 ±0.5 

38 ±3 

16.5 ±0.2 

8.7 ± 2 

4.4 

140 

35 

G138-52 

18.5 ±0.5 

80 ±6 

18.6 ±0.2 

80 ±6 

1.0 

52 

38 

G198+32 

20.0 ±0.5 

53 ±9 

18.4 ±0.2 

16 ±2 

3.2 

60 

12 

G198+32 E 

20.0 ±0.5 

80 ±29 

20.0 ±0.2 

138 ± 19 

0.6 

7 

7 

G221+35 

19.1 ±0.4 

97 ±5 

18.3 ±0.2 

48 ±3 

2.0 

32 

22 

G236+39 

19.5 ±0.6 

74 ± 17 

18.4 ±0.2 

42 ±2 

1.8 

45 

30 

G236+39 N 

19.5 ±0.6 

51 ±5 

19.0 ±0.5 

54 ±6 

1.0 

28 

22 

G254+63 

20.5 ±0.8 

121 ± 14 

18.9 ±0.2 

33 ±2 

3.6 

84 

27 


[G/D\ lo , and [G/D] diff are gas-to-dust mass ratio inferred from the slope of dust optical depth versus H I column 
density for the clouds and for the regions just outside the clouds, respectively. For reference the average interstellar 
medium value is 124 (Li & Draine|[2001). The ‘dark’ factor, /dark is defined in Eq. [ t] 








3 4 5 6 7 8 9 

N(HI) from Arecibo [10 20 cm' 2 ] 

Fig. 2.— (Top panel) An image covering 8.2° x 7.5° centered on J2000 coordinates 23:11:51 +21:11 
(galactic coordinates 93.9, -36.0) with celestial North upward. The color composite was constructed 
with the gas column density (inferred in the optically-thin limit) from the Arecibo 21-crn line in 
red+blue and the dust column density inferred from the Planck optical depth in green. (Bottom 
panel) Hess diagram (densitized scatter plot) of the dust versus gas column density using the pixel 
values from the image shown in the top panel. The density of points in the scatter plot is shown 
in a continuous color table from red (high-density) to blue (low-density). 
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3 4 5 6 7 8 9 

N(HI) from Arecibo [10 20 cm' 2 ] 

Fig. 3.— (Top panel) An image covering 8.2° x 7.4° centered on J2000 coordinates 23:44:14 +08:17 
(galactic coordinates 95.9, -51.0) with celestial North upward. The color composite was constructed 
with the gas column density (inferred in the optically-thin limit) from the Arecibo 21-crn line in 
red+blue and the dust column density inferred from the Planck optical depth in green. (Bottom 
panel) Hess diagram (densitized scatter plot) of the dust versus gas column density using the pixel 
values from the image shown in the top panel. The density of points in the scatter plot is shown 
in a continuous color table from red (high-density) to blue (low-density). 



3.5 



0.0 _I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_ 

3 4 5 6 7 8 

N(HI) from Arecibo [10 20 cm' 2 ] 

Fig. 4.— (Top panel) An image covering 7.1° x 6.4° centered on J2000 coordinates 23:46:35 
+22:17 (galactic coordinates 103.7, -38.2) with celestial North upward. The color composite was 
constructed with the gas column density (inferred in the optically-thin limit) from the Arecibo 21- 
cm line in red+blue and the dust column density inferred from the Planck optical depth in green. 
(Bottom panel) Hess diagram (densitized scatter plot) of the dust versus gas column density using 
the pixel values from the image shown in the top panel. The density of points in the scatter plot 
is shown in a continuous color table from red (high-density) to blue (low-density). 
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N(HI) from Arecibo [ 10 20 cm' 2 ] 

Fig. 5.— (Top panel) An image covering 8.0° x 7.4° centered on J2000 coordinates 00:15:09 
+08:53 (galactic coordinates 108.0, -52.9) with celestial North upward. The color composite was 
constructed with the gas column density (inferred in the optically-thin limit) from the Arecibo 21- 
cm line in red+blue and the dust column density inferred from the Planck optical depth in green. 
(Bottom panel) Hess diagram (densitized scatter plot) of the dust versus gas column density using 
the pixel values from the image shown in the top panel. The density of points in the scatter plot 
is shown in a continuous color table from red (high-density) to blue (low-density). 
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3 4 5 6 7 8 9 

N(HI) from Arecibo [10 20 cm' 2 ] 


Fig. 6.— (Top panel) An image covering 7.1° x 6.4° centered on J2000 coordinates 01:27:51 
+09:40 (galactic coordinates 137.7, -52.2) with celestial North upward. The color composite was 
constructed with the gas column density (inferred in the optically-thin limit) from the Arecibo 21- 
cm line in red+blue and the dust column density inferred from the Planck optical depth in green. 
(Bottom panel) Hess diagram (densitized scatter plot) of the dust versus gas column density using 
the pixel values from the image shown in the top panel. The density of points in the scatter plot 
is shown in a continuous color table from red (high-density) to blue (low-density). 
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N(HI) from Arecibo [10 20 cm' 2 ] 

Fig. 7.— (Top panel) An image covering 8.3° x 5.8° centered on J2000 coordinates 08:29:24 
+26:29 (galactic coordinates 197.3, +32.2) with celestial North upward. The color composite was 
constructed with the gas column density (inferred in the optically-thin limit) from the Arecibo 21- 
cm line in red+blue and the dust column density inferred from the Planck optical depth in green. 
(Bottom panel) Hess diagram (densitized scatter plot) of the dust versus gas column density using 
the pixel values from the image shown in the top panel. The density of points in the scatter plot 
is shown in a continuous color table from red (high-density) to blue (low-density). 
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2 3 4 5 6 V 0 

N(HI) from Arecibo [10 20 cm' 2 ] 

Fig. 8.— (Top panel) An image covering 8.1° x 5.9° centered on J2000 coordinates 09:10:15 
+09:01 (galactic coordinates 221.1, +34.9) with celestial North upward. The color composite was 
constructed with the gas column density (inferred in the optically-thin limit) from the Arecibo 21- 
cm line in red+blue and the dust column density inferred from the Planck optical depth in green. 
(Bottom panel) Hess diagram (densitized scatter plot) of the dust versus gas column density using 
the pixel values from the image shown in the top panel. The density of points in the scatter plot 
is shown in a continuous color table from red (high-density) to blue (low-density). 
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N(HI) from Arecibo [10 20 cm' 2 ] 


Fig. 9.— (Top panel) An image covering 7.7° x 4.5° centered on J2000 coordinates 09:47:04 
+02:47 (galactic coordinates 233.7, +39.7) with celestial North upward. The color composite was 
constructed with the gas column density (inferred in the optically-thin limit) from the Arecibo 21- 
cm line in red+blue and the dust column density inferred from the Planck optical depth in green. 
(Bottom panel) Hess diagram (densitized scatter plot) of the dust versus gas column density using 
the pixel values from the image shown in the top panel. The density of points in the scatter plot 
is shown in a continuous color table from red (high-density) to blue (low-density). 



16 - 



0 2 4 6 8 10 


N(HI) from Arecibo [10 20 cm' 2 ] 

Fig. 10.— (Top panel) An image covering 8.0° x 7.4° centered on J2000 coordinates 11:30:40 
+07:51 (galactic coordinates 254.4, +62.8) with celestial North upward. The color composite was 
constructed with the gas column density (inferred in the optically-thin limit) from the Arecibo 21- 
crn line in red+blue and the dust column density inferred from the Planck optical depth in green. 
(Bottom panel) Hess diagram (densitized scatter plot) of the dust versus gas column density using 
the pixel values from the image shown in the top panel. The density of points in the scatter plot 
is shown in a continuous color table from red (high-density) to blue (low-density). 
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Several peaks within the clouds are listed in the Planck Galactic Cold Core Catalog (PGCC), 
which contains all compact cores from the Planck survey with dust temperature significantly colder 
than surroundings (Planck Collaboration et al. 2015). Table [3] lists the cataloged cold cores within 
our sample clouds together with the clump dust temperature and emissivity index. The weighted 
mean dust temperature is (T c i ump ) = 14.85 ± 0.07 K, and the weighted mean emissivity index is 
(/3) = 1.81 ± 0.31. While the (3 values are comparable to those we estimated from entire clouds, 
the temperatures in the PGCC are colder. This measures the decrease of dust temperature and 
increase of /3 that continues from the diffuse medium (average Td = 20.8 K , /3 = 1.55 over the 
faintest 10% of the sky (Planck Collaboration et al. 2014])) to the isolated clouds ((Td) = 18.1 ±0.1, 
((3) = 1.66 from Table [2]) down to the smaller scale of the cold cores. For some of the cold cores, 
there are distance estimates in the PGCC based on reddening of SDSS stars. For G94-36, the 
distances are in the range 203-281 pc (and one at 421 pc); for G198±32 the distance estimates are 
290-360 pc; and for G103-39 the distances estimate is 241 pc. The distance estimates should be 
treated with caution because if the low extinctions and uncertainties in stellar distances, but the 
suggest at least some of the clouds are further than the fiducial 100 pc we would normally assume 
for high-latitude interstellar clouds within the scale height of cold atomic gas. 


3. Results 


3.1. Properties of the clouds 


For each cloud, the dust temperature was measured by Planck and IRAS using Equation [T] 
In no case was one of the clouds of our sample warmer than its surroundings, indicating that none 
of them are heated internally by stars. Table [2] lists the dust temperatures, Tdiff, in the diffuse 
medium around the clouds as well as the dust temperatures, T c ioud) i n the brightest parts of the 
clouds. The dust temperatures in the clouds are typically 1 K lower than the surroundings, with 
temperature decrements into the clouds spanning 0 to 1.6 K. 


Of critical importance for this project is the manner in which the dust optical depth was 
derived from the observed intensities. The dust column density is derived from multi-frequency 
measurements that cover the spectral energy distribution from 100 to 850 /irn, spanning both 
sides of the peak and including sufficiently long wavelengths to measure the dust temperature and 
wavelength-dependent emissivity. The emissivity index, (3, was determined using maps degraded in 
resolution to 30 / , in order to improve sensitivity, while the optical depths were then calculated at 
the full 5 ; resolution. This will lead to systematic errors in locations where /3 is spatially variable. In 
fact, the clouds in this study tend to be local maxima in /3. If (3 were measured at full resolution, the 
maxima would likely become more pronounced. Because the dust temperature and (3 are inversely 
correlated (Juvela et al. 2013), the dust temperatures in the clouds would likely decrease if (3 were 
measured at full resolution. For these reasons, we expect the cloud temperatures are likely to be 
even lower than the surroundings. The amount of dust required to explain the observed emission 
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Table 3: Cold Cores within the Isolated Clouds 


Cloud 

PGCC a 

-^clump 

/^clump 



(K) 



G104-39 G105.76-38.36 

G106.06-38.39 
G105.07-38.06 
G103.35-38.93 
G103.38-39.33 
G108-53 G108.74-52.67 

G109.26-52.52 
G198+32 G198.25+32.30 

G197.22+31.99 
G197.37+31.73 
G197.91+31.99 
G198.40+32.26 
G196.55+31.96 
G221+35 G220.54+34.70 

G254+63 G253.01+61.23 

G94-36 G94.19-37.37 

G91.71-34.32 
G94.36-36.31 
G92.07-36.71 
G94.28-36.61 
G92.34-34.41 
G94.04-34.44 
G92.59-35.55 
G92.06-34.59 
G94.68-37.21 
G92.38-37.43 
G93.95-34.17 
G92.23-35.70 
G96-51 G96.05-50.34 

G95.72-49.63 

G236+39 G235.60+38.28 

G235.67+38.00 


14.8 

± 

1.7 

1.98 

± 

0.38 

14.7 

± 

2.5 

1.79 

± 

0.53 

15.3 

± 

1.5 

1.96 

± 

0.27 

14.0 

± 

2.4 

1.93 

± 

0.62 

14.9 

± 

1.9 

2.04 

± 

0.39 

15.5 

± 

1.5 

1.78 

± 

0.29 

13.5 

± 

3.6 

2.19 

± 

1.11 

17.1 

± 

3.1 

1.22 

± 

0.61 

15.6 

± 

3.4 

1.40 

± 

0.80 

16.1 

± 

3.6 

1.55 

± 

0.79 

14.9 

± 

3.2 

1.65 

± 

0.74 

18.3 

± 

2.5 

1.07 

± 

0.38 

18.5 

± 

3.5 

1.29 

± 

0.64 

13.7 

± 

1.9 

2.34 

± 

0.44 

14.7 

± 

3.5 

1.44 

± 

0.86 

13.5 

± 

1.3 

1.83 

± 

0.28 

14.0 

± 

1.4 

2.06 

± 

0.29 

14.1 

± 

1.6 

2.12 

± 

0.40 

12.6 

± 

1.2 

2.27 

± 

0.34 

15.6 

± 

2.4 

1.83 

± 

0.45 

13.4 

± 

1.8 

2.08 

± 

0.49 

15.2 

± 

2.5 

1.80 

± 

0.52 

15.3 

± 

2.5 

1.85 

± 

0.59 

13.1 

± 

3.5 

2.33 

± 

1.06 

16.9 

± 

2.8 

1.61 

± 

0.52 

16.5 

± 

2.4 

1.73 

± 

0.42 

15.0 

± 

3.0 

1.96 

± 

0.74 

16.2 

± 

2.9 

1.71 

± 

0.52 

17.7 

± 

1.8 

1.39 

± 

0.23 

17.5 

± 

2.5 

1.67 

± 

0.46 

14.6 

± 

1.9 

2.11 

± 

0.46 

16.5 

± 

2.2 

1.91 

± 

0.42 


Designation within the Planck Galactic Cold Clump Catalog (Planck Collaboration et al.||2015|) 
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increases if the dust temperature is lower. Therefore, for the same observed emission from the 
clouds, we expect a full-resolution determination of /3 and T c i ou d would yield dust optical depths 
higher than we report. The amount by which f3 might increase is not known but is likely no more 
than 0.2 (based on the observed range of (3 on large scales and not expecting drastic spectral 
variations on angular scales between 5 7 and 30'), with implied systematic uncertainty of the cloud 
optical depths at the 30% level. We emphasize this point because we are measuring dust excess 
relative to gas in clouds, and a full-resolution measurement would tend to make the dust excesses 
even more pronounced than we measure. 


To give some context for the clouds, Table [2] lists the mass of each cloud, measured using the 
maps of t(353 GHz) and the standard gas-to-dust mass ratio of 124 (Li & Draine 2001), and the H I 
mass, measured using the 21-cm line integral and assumed optically thin. Both mass determinations 
were made in 80 7 radius regions (grown to 100 7 for the largest cloud) centered on the infrared peak, 
with the unrelated emission in front of and behind the cloud estimated in an annulus surrounding 
it. We assume the clouds are at a distance 100 pc, typical of the scale height for cold atomic gas 
(Dickey & Lockman 1990). The clouds are typically 100 Mq. 


3.2. Comparison between dust and gas 


Figures 10 through [3] summarize the observations. The top panel of each Figure shows a 
color composite image combining the Arecibo and Planck data. The fields are filled with emission 
from the diffuse interstellar medium, with relatively bright features, which we noticed from the 
all-sky images as relatively isolated clouds, near the center. The lower panels show pixel-by-pixel 
comparisons of the images in the form of Hess diagrams, wherein we count the number of image 
pixels that fall within each of 300 x 300 cells in the 1V(HI), r(353 GHz) plane. For the lowest-column 
densities, which pervade much of the area of the images, and hence most of the points in the Hess 
diagrams, there is a fairly tight correlation between dust and gas. This correlation is evident by the 
relatively constant color of the spatial images, and the high density of points forming a red ridge 
in the Hess diagrams. 


For the distinct cloud-like features, the dust and gas trace similar structures but with distinctly 
different relative amounts. For example, in Figure |To| the ellipses labeled A, B, C show the spatial 
location (top panel) and the dust-gas correlation (bottom panel) of three features. The coherence 
of the features in the image and in the Hess diagram proves that the features have different relative 
amounts of dust and atomic gas. It is evident that within the cloud-like features, there is relatively 
more dust per atomic gas than there is in the more diffuse medium. The individual clouds in the 
top panel correspond to ‘tracks’ in the bottom panel due to their different dust-gas proportions. In 
Figure [5j the circles labeled A and B tell similar stories. Because the ‘green’ emission is associated 
with the structure of the H I cloud, it is almost certainly all interstellar. But because galaxies could 
also cause infrared excess (they gas being at velocities outside the 21-cm survey), we investigated 
the relatively compact features in the ellipse labeled B individually. The easternmost ‘green’ feature 
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in ellipse B is at 00 h 30 m 58.9 s +10°46 / 30 // , and there are no corresponding sources in the SIMBAD 
database. The peak just to its west is at 00 h 29 m 38.1 s +10°47 , 36 ,/ and has a corresponding IRAS 
point source catalog entry at 100 //rri only with a flux density of 3 Jy. Inspection shows the source 
is extended and there is no optical counterpart; the most likely explanation is that the IRAS and 
Planck emission are both interstellar. 


3.3. Measuring the gas/dust ratio 


The Figures of the clouds show that there are compact (‘green’) regions where there is material 
being traced by the dust maps but not by the 21-cm maps. We can quantify the ‘excess’ emission by 
scaling and subtracting the 21-cm map. The measurement procedure is best summarized with an 
example. Figure [TT] shows slides through two components of the same cloud, G236+39. Note how 
the northern component has a simple, linear trend of dust versus gas. The material along this slide 
is most likely all atomic, and the 21-cm line is optically thin, so that both r(353 GHz) from Planck 
and IV(H I) from Arecibo are linear tracers of column density of the same material. The slope of 
this line is the nominal dust optical depth per unit gas column density (r/IV(H)). In contrast to 
the northern, purely atomic region, the cut through the southern component of the cloud shows a 
marked infrared excess at higher column densities. The linear fit from the northern component is 
still valid for the outermost portions of the southern component, forming a lower envelope at the 
higher column densities. Now examine the spatial profiles of these clouds, in the right-hand panels 
of Figure 11 While the dust and gas are a very good match for the northern component, leaving 
little residual after subtracting the gas from the scaled dust map, the southern component shows a 
significant positive excess, which is shown as the thin solid curve and is defined as 

r(353 GHz) 


= 


- IV(H I). 


(2) 


(r/N( H)} 

In words, the infrared excess column density is the dust optical depth, scaled by the slope of the 
dust/gas correlation, then subtracting the atomic gas. In Figure 11, the excess can be directly 
measured from the right-hand panels by taking the positive deviation of the dust optical depth 
above the straight line fitted at low column densities, then dividing that deviation by the slope 
of the line fitted at low column densities. The dust opacity per unit gas column density and the 
infrared excess column density can now be derived qualitatively, to interpret the color images of 
the clouds. 


For each cloud, we identified the location of the infrared excess peak in the Planck r(353 GHz) 
map, and generated the H I 21-cm line spectrum from the Arecibo GALFA survey (subtracting a 
nearby off-cloud position to remove the large-scale HI). In all cases, the velocity component in the 
H I spectrum associated with the dust peak could be readily identified. The central velocity and 
the full-width-at-half-maximum (FWHM) of these 21-cm line components is listed in Table [I] We 
then generated an H I map covering the velocity channels centered on that line profile. These H I 
maps in one velocity range had a much better correspondence to the Planck images, for the cloud 
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centered on the infrared excess peak. We then generated a one-dimensional profile through each 
cloud, using a sliding box of width 1° for each 1' pixel. The dust and velocity-selected H I profiles 
showed peaks at the appropriate locations, and the slope of the gas-to-dust trend for these profiles 
was computed with a linear fit. 


The dust opacity per unit gas column density is directly related to the gas/dust mass ratio. 
The dust optical depth 

t u = K l/ N(R)fj,mu (3) 

where mu is the mass of a hydrogen atom, ji is the gas mass per unit H-atom taking into account 
other elements (/r = 1.4), and -/V(H) is the hydrogen column density. Note that this is the total 
column density of H nuclei, including, if they exist and have associated dust, molecular and ionized 
gas. The mass opacity coefficient, depends only on the material properties and size distribution 


of the dust grains. Using k(353 GHz; 850 /xm) = 0.42 cm 2 g 1 from Li Sz Draine (2001), we can 
directly convert the observed dust opacity per unit gas column density into a gas/dust mass ratio: 


[G/D] = 9.7 x 10" 25 / 


r(353GHz) 
iV(H I) 


( 4 ) 


where the H I column density is expressed in cm~ 2 . The gas/dust mass ratios for the clouds are 
listed in Table [l] as [G/-D] cloud . To measure the gas/dust ratio for regions surrounding the clouds, 
we took the total H I column density, minus the component centered on the cloud, and measured 
the gas/dust slope across the 64 deg~ 2 region around each cloud, using a robust fit that excludes 
positive outliers. These values are listed in Table [2] as [G/D] diS . Please keep in mind that the 
[G/D\ values assume all the gas is atomic. We discuss possible interpretations below, but a low 
value of [G/D] generally means means there is a deficiency in the inventory of gas. 


It is evident that the gas/dust ratios are significantly higher in the diffuse medium than in 
the clouds. In the diffuse interstellar medium, the wide-held observations by COBE used in the 
Li & Draine (2001) model yield a gas/dust mass ratio of 124. This value is an “upper envelope” 


for the regions in the present Arecibo /Planck study that was centered on infrared excess peaks. 
The cloud cores all have gas/-to-dust mass ratios significantly lower than the diffuse interstellar 
medium. Many of the regions surrounding the clouds have lower gas/dust ratio than the diffuse 
ISM, suggesting that the same effect that causes the infrared excess peaks spreads over significant 
fractions of the regions studied. 
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Fig. 11.— Comparison of dust and gas in two horizontal slices through G236+39. The upper 
panels are for a cut through the northern portion, G236+39 N, at constant declination (J2000) 
+02°04 / , while the lower panels are for a cut through the southern portion, G236+39 S, at constant 
declination +00°32 / . (left) The H I could density (dashed) and dust optical depth (thick solid) as 
a function of offset in right ascension along the slices. Notice how for G236+39 N (upper panel), 
the gas and dust follow each other closely. In contrast, notice how for G236+39 S (lower panel), 
there is a much sharper peak of the dust optical depth than there is of the gas column density. 
(right) Scatter plots of the dust optical depth versus H I column density using the data from the 
left-hand panels. For G236+39N (upper panel), where the dust and gas are well correlated, a simple 
linear fit (thin solid line) fits the data showing the dust and atomic gas are well mixed and linearly 
correlated. A red line shows a median B-spline fit, while the green lines surrounding it are the 
25% and 75% percentile splines. For G236+39 S (lower panel), we over plot the same lines from 
G236+39 N, without scaling; this shows that at low column densities the gas-to-dust ratio is the 
same as in G236+39 N, while at higher column densities there are very large positive deviations 
indicating excess dust. 
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4. Discussion 
4.1. Trends of Dust/Gas 



Fig. 12.— The gas/dust mass ratio, measured from the optical depth at 353 GHz per unit H I 
column density, versus the dust temperature. The cloud peaks are shown as asterisks (*), and the 
diffuse regions around the clouds are shown as diamonds (o). 


The gas/dust ratio varies systematically with the temperature of the dust. Figure 12 shows 
the gas/dust ratio versus dust temperature, combining the observations toward the cloud peaks 
and the diffuse medium near the clouds. The diffuse medium dust temperatures were obtained 
from modes of the histograms of the 8° x 8° regions surrounding the cloud. The clouds have lower 
gas/dust ratio compared to the interstellar gas-to-dust ratio of 124 (Li Sz Draine 2001); the cloud 
peaks tend to be lower than the regions surrounding them diffuse medium; and the colder clouds 
and diffuse regions have lower gas/dust than the warmer ones. The dynamic range of the trend is 
an order of magnitude. 
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The apparently excess dust is not an artifact of converting the observed infrared intensities 
into the dust optical depth. The dynamic range of the dust temperature variations is only 17-20 
K, over which the predicted emission at 353 GHz changes by only 28%, which is far smaller than 
the change in gas/dust. More likely, the dust temperature in the cloud cores is overestimated 
because of temperature gradients. The dust temperatures in the cloud cores are lower than that of 
their surroundings. If there is a radial temperature gradient into the core, the central temperature 
is lower than the apparent temperature of dust as seen by us, because both the core and the 
warmer dust surrounding it are present on the same line of sight. In this case, the optical depth 
is underestimated and the actual gas/dust in the clouds is actually even lower than inferred from 
the observed line-of-sight averages. The decreased gas/dust ratio in the clouds is not due to the 
conversion from infrared intensity to dust column density, and we turn to physical explanations for 
the effect. 


Hypotheses for the observed trend are (1) the gas content neglects to include ‘dark’ molecular 
gas (14.2); (2) the amount of atomic gas in the clouds is underestimated because the 21-cm line is 
and (3) the dust properties change such that the optical depth per unit dust 


optically thick 

mass is higher in colder clouds( ^4.4[ ). We discuss each of these hypotheses in turn in the following 
subsections. 


4.2. Interpreting as ‘Dark Gas’ 


If we interpret the dust emission as a tracer of the total gas column density, using the correlation 
in the diffuse medium, then inverting Eq. [4j 

iV(H) dust = 1.03 x 10 24 r(353GHz)(G/D) difr , (5) 

where (G/D) is the gas-to-dust ratio measured measured in a reference region near the cloud or 
taken from an interstellar average ({G/D)ism = 124). If the low values of gas/dust mass ratio are 
due to dust associated with some “dark” gas, then the column density of that gas 

iV(H) dai . k = IV(H) dust - IV(H I). (6) 


Note that this column density is the equivalent number of nucleons, as for atomic gas; if the gas 
were molecular, then lV dark = 2JV(H2). The factor by which the total gas column density exceeds 
the atomic gas is then 


IV(H) tot 

J dark — 


[G/D] 


diff 


(7) 


iV(HI) [G/D] cloud ’ 

measures the amount of ‘dark’ gas with a unitless quantity that is directly (and linearly) related 
to observables. (Note that the dust properties needed for these calculations, k, cancel out in Eq. [7] 
as long as they remain constant from cloud to more diffuse ISM.) If / dar k = 1 then the cloud has 
the same gas/dust as the surrounding medium (and no ‘dark’ gas is required), while larger values 
indicate the ratio of total column density to that of atomic gas. The values of / dark are summarized 
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in Table [2j For all clouds but one, the inferred amount of ‘dark’ gas is larger than that of atomic 
gas, with total masses exceeding atomic gas by factors 2-4. 0 The fraction of the total amount of 
gas that is ‘dark’ is 1 — * k , with typical values 50% to 75%. In smaller patches comparable to 

the 5 7 beam size, the relative amount is significantly higher, reaching up to 6, meaning dark gas 
comprises 83% of the total. The ratio of ‘dark’ mass integrated over the clouds can be estimated 
similarly, though without taking into account the gas-to-dust ratio of the environment, using the 
cloud mass determined from the dust and H I maps; values of M c i ou d/AfH I are similar to fdark in 
the range 1-5. 

The most likely candidate for ‘dark’ gas is H 2 . Molecular hydrogen forms readily in interstellar 
conditions and would dominate the ISM if it were not dissociated by ultraviolet photons from 
the diffuse interstellar radiation field. In this hypothesis, IV(H 2 ) = ^IVdark measures the column 
density of molecular gas. This hypothesis can be tested in many ways: comparison to existing 
tracers of molecular gas, comparison to results of absorption-line observations of molecular gas 
along essentially random lines of sight (toward ultraviolet sources), and searches for new tracers of 
molecular gas. The last idea is the subject of follow-up observations we are pursuing, using OH 
absorption and [C II] emission. We address each of the other two tests now. 


4-2.1. Comparison to CO 


For high galactic latitude regions, there is no sensitive all-sky survey of the best mapping tracer 
of molecular gas, which is the millimeter-wave CO(l-O) line. However it was found that the Planck 
100 GHz band contained significant emission from the CO(1-0) line that could be separated from 
the continuum emission using a multi-frequency analysis (Planck Collaboration XIII 2014). The 
Type 3 CO map from that paper has a resolution of 5.5 7 and rms fluctuations of the (mostly empty) 
map near one of our selected clouds are 0.5 K km s” 1 at this angular scale. For molecular clouds 
in the galactic plane, the typical ratio X(CO) = iV(H 2 )/W(CO) = 2 x 10 20 cm” 2 K” 1 km” 1 s, 
so the Planck CO map would be sensitive to iV(H 2 ) ~ 10 20 cm” 2 if that same ratio applied. For 
comparison, the map of IVdark has an rms at 6 7 angular scale of 0.4 x 10 20 cm” 2 , which would 
correspond to molecular column densities of 0.2 x 10 20 cm” 2 . Thus the 5 a clouds in the Planck 
maps of IVdark should appear in the CO map, if they have the same scaling as galactic plane 
molecular clouds. 


The Planck CO map for some clouds shows a good correspondence with the infrared excess, 
in particular G94-36, G104-39, and G108-53; others show a CO detection over a limited portion 
of the infrared excess region, such as G96-51, G138-52, and G236+39; while the rest show no CO 


1 In many cases, one should consider fdark calculated in this manner as lower limits, because the ‘dark’ gas may 

extend far from the cloud and into the region where [G/Z?]diff is calculated. We use the conservative values in the 

remainder of the paper, but note for now that if instead we used the diffuse ISM gas/dust ratio value, (G/D) ism as 

the reference, the values of /dark would be typically 3-8. 
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Fig. 13.— The ‘dark’ nucleon column density versus the CO line integral for the cloud G108- 
53. They color image is the Hess diagram showing the density of pixels at each value of A^ark 
and VF(CO), generated with a gaussian kernel with width 0.5 in each axis (in the units shown). 
The straight dashed line is the slope that would be expected based on the galactic plane value of 
Xqo = 2.1, while the dotted line is a rough lower bound at Xd ar k/kF(CO) ~ 1, which corresponds 
to X co = 0.5. 


detection. Figure [13] shows a pixel-by-pixel comparison of the W( CO) and N^k for G108-53. The 
slope expected for the X-factor of the galactic plane molecular clouds forms a rough upper bound 
of the points, while a lower bound is formed by a line with a 4 times lower effected X-factor. The 
CO and ‘dark’ gas maps agree in detail for this cloud, and after subtracting the CO image scaled 
by the dashed line, the prominent infrared cloud that is green in Figure [5] is removed cleanly. To 


quantify the correlation between CO and infrared excess, we made Hess diagrams like Figure 13 


for each cloud and estimated the slope |Xd ar k/bF(CO). (The factor of \ is needed because we 
calculate the dark gas column density in units of H nuclei, while the X-factor is for H 2 molecules.) 
For many clouds, only lower limits were possible due to non-detection of CO. Then we identified 
the peaks in the infrared excess and the CO map (when present in the latter) and measured the 
ratio directly on the peaks. Again, some ratios are lower limits due to lack of detected CO, and we 
use 95% confidence (2 a) limits. Table [4] summarizes the results. 
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Table 4: “X-factor” for Infrared excess clouds 


Name 

(^darkWCO))* 

^dark/^(CO) b 

G94-36 

0.45 ±0.05 

0.34 ±0.05 

G96-51 

2.3 ±0.4 

2.8 ±0.3 

G104-39 

0.60 ±0.1 

0.60 ± 0.1 

G108-53 

0.75 ±0.1 

0.65 ±0.1 

G138-52 

1.1 ±0.3 

0.90 ± 0.1 

G198+32 

> 4 

> 5 

G221+35 

> 3 

> 4 

G236+39 

1.2 ±0.2 

1.4 ±0.1 

G254+63 

> 10 

> 20 


“half the slope of the infrared excess column density, 7Vd ar k, versus CO line integral, WTCO), from the Hess diagram 
^half the ratio of IVdark over CO line integral W(CO) at the cloud peak 


For the brightest infrared-excess clouds with clearly associated CO peaks, the values of X (CO) 
we find are about a factor of 2 lower than the galactic plane value. For G94-36 there is a linear 
correlation between IVdark and VF(CO), yielding an accurate measure of X(CO). Figure 13 shows 
for G108-53 a range of apparent slopes between IVdark and W(CO), possible curving from a high 
slope in faint regions to a lower slope in bright ones (though this cannot be quantified accurately 
due to the low signal levels). For other clouds, the lower limits on X(CO) exclude the galactic 
plane value (cf. Bolatto et al. 2013) as well as the value for the brighter infrared excess clouds 
from this paper. The high X(CO) for those clouds does not rule out the presence of H 2 with faint 
CO. The linear scaling between CO brightness and H 2 column density only applies to molecular 
clouds that are gravitationally bound (Bolatto et al. 2013). If the CO line is optically thin or sub 
thermally excited, X(CO) can be high even if CO is abundant, so we cannot rely upon the CO(l-O) 
emission as a tracer in diffuse regions. 


4-2.2. Trend of dark gas column density with dust temperature 

The inferred amount of ‘dark gas’ depends very cleanly on the difference in dust temperature 
between the cloud and its surroundings. Figure 14 shows inferred amount of ‘dark gas’ increases 
monotonically with the dust temperature difference between the cloud interior and surroundings. 
This trend has a relatively straightforward interpretation. Within cold cloud interiors, molecules 
form and the H I is no longer tracing the total column density. In the context of the ‘dark gas’ 
hypothesis, the cold regions are where the dust is shielded from the radiation field, and those 
are the same regions where H 2 self-shields. As the amount of self-shielding increases the dust 
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^cloud ^ diff 1^1 


Fig. 14.— The ratio of total column density to atomic column density, /dark, for each cloud versus 
the dust temperature difference between the interior of the cloud and its surroundings. 


temperature decreases, and the molecular hydrogen formation rate surpasses the dissociation rate 
allowing hydrogen to become largely molecular. The regions with little self-shielding have dust 
temperature as warm as the widespread ISM and the gas remains photo dissociated. 


The trend between dark gas and dust temperature can be reproduced by a simple model. The 
abundance of H 2 depends on the abundance of H I from which it forms, as well as the dissociation 
rate; for a clod with uniform density 100 cm -3 and gas temperature 80 K, 


21V(H 2 ) 
IV(H I) 


0.8 X - 2 iV(H I), 


( 8 ) 


where x is the amount of radiation relative to the diffuse interstellar radiation field (Reach et al. 


1994 


eq. 3). The energy emitted (or absorbed) by dust depends on T 4+/3 , so the dust temperature 


difference is related to the radiation field by 


T = T 0 x 4 +e 


( 9 ) 


where To is the dust temperature at the surface of the cloud. The attenuation of the radiation field 
is roughly approximated as 


X = e 


2A V /1.086 


( 10 ) 
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where the visible-light extinction is 

A v = [JV(H I) + 2JV(H 2 )] /1.9 x 10 21 cm' 2 


( 11 ) 


(Savage Sz Mathis 1979). The factor of 2 in the equation for % is to account for the absorption of 


dust-heating photons being at ultraviolet rather than visible wavelengths. Note that the gas-to- 
dust ratio in Equation [Tl] is consistent with the assumed value for the diffuse interstellar medium 
by construction because it was part of the dust model from which obtained k. The prediction of 


this simple model is shown in Figure 14 and agrees quite well with the observed data. Thus the 


idea of the infrared excess being due to dust associated with ‘dark’ molecular gas can match the 
observations quantitatively. 


4.3. Interpreting as optically thick H I 


We have assumed throughout this paper that the integral over the 21-cm line brightness traces 
the column density of the atomic gas linearly. In fact, the 21-cm line is prone to optical depth 
effects over the column densities of interest. The under prediction of the H I column density by the 
optically thick 21-cm line could account for the apparent excess dust column density in interstellar 
clouds, as was recently suggested by e.g. Fukui et al. (2015). The optical depth 


t(v) = N(v)/CT S 


( 12 ) 


where v is the velocity, N(v) is the column density per unit velocity at velocity v, C = 1.823 X 
10 18 cm -2 KT 1 km -1 s, and the line brightness 


T B {v) = T s [ l-e- T M], 


(13) 


where T s is the excitation (‘spin’) temperature (Kulkarni & Heiles 1988). The column density in 
the optically thin limit 

N* H = C [ T B {v)dv 


while the actual column density 


Nh = J N{y)dV = CT S j r{v)dv. 


To calculate the total column density, one must correct the column density in the optically thin 
limit by a factor 

N h / r \ _ / T B {v)^^dv ^ 


( 15 ) 


N* h \ 1 — e~ T / f T B (y)dv 

For a single gaussian component with peak optical depth To, we find a good fit 

0.75 


Nh 
N 


H 


TO 


1 - e _T ° 
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over a wide range of optical depths (0.01 < tq < 20) and independent of spin temperature. 


An advantage of studying isolated clouds at high galactic latitude is that the observations are 
dominated by a single velocity component of 21-cm emission and a single spatial feature in both 
21-crn and infrared emission. If the clouds we are observing are entirely atomic gas, and N& us t, as 
in the previous section, measures the total column density, then the ratio 


Adust 

N 


— /dark 


(16) 


H 


which can be combined with Equation [15] to solve for the optical depth of the 21-cm line. Taking 
the observed /dark for our clouds from Table [2] the implied 21-cm line optical depths of the clouds 
are of order unity, with typical values tq ~ 3. If the 21-cm lines are indeed this optically thick for 
the majority of the clouds, then observations of radio sources will show deep absorption features. 
The results of the Millennium H I survey (Heiles & Troland 2003) indicate that optical depths of the 
21-cm line are typically small, but a pencil-beam survey cannot be used to address the possibility 
of regions of high 21-cm optical depth that coincide with where we observe high gas/dust ratios. 
The possibility that the apparently low gas/dust ratio in isolated clouds could be due to the H I 
being optically thick can be followed up with future observations of radio sources situated behind 
the clouds, and will be the subject of future observational follow-up. 


In a recent paper, Fukui et al. (2015) suggest that the excess infrared emission often referred 


to as ‘dark gas’ is optically thick H I, which is the same hypothesis we explore in this section of 


this paper. Fukui et al. (2015) measure the temperature variation of dust using Planck and find 


it is correlated with the apparent infrared excess, then use this empirical correlation to determine 
the total column density under the assumption the extra infrared emission is from optically thick 
H I. Shielding of the radiation field from the interior of clouds is a straightforward explanation 
that must be contributing to the apparent trends. The simple model we presented in the previous 
section (Eq. [8H9]) and illustrated in Figure 14 seem to match the observations for the clouds in 


our sample well. The lack of widespread optically thick H I from 21-cm absorption studies further 
argues against significant an amount cold H I adequate to explain the range of gas/dust variations 


that are observed. On a technical note, our equations above (Eq. 14 15) are similar to theirs; 


however, they authors did not take into account the line profile and defined a single optical depth 
for all velocities along each line of sight. This fictitious optical depth , t*, can be solved from the 
line profile and Equation |14[ 


1 — e~ 


1 — e~ 


Fukui et al. (2015) use the fictitious r* together with a fictitious line-width AV (the ratio of the line 
integratal to the peak brightness temperature), and a single spin temperature T s . as intermediaries 
to convert the observed line brightness into a total column density. None of r*, AV, or T s apply on 
typical lines of sight where the WNM contains about half the column density and has a distinctly 
different profile than the CNM. The Millennium Survey results showed that randomly selected 
lines of sight toward extragalactic sources have multiple velocity components of varying widths and 
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temperatures (Heiles Sz Troland 2003). For the simplest case of a single CNM line on the line of 


sight, like for the clouds in the present paper, the concepts are appropriate, and the more detailed 


solution finds that Fukui et al. (2015) would systematically overestimate the true column density 


by an amount that varies from 0 at low optical depth, to 15% at to = 1, to 50% at to = 5, where 
to is the peak optical depth of the line profile. 

Our results being produced solely by optically thick 21-cm line is highly unlikely because it 
requires typical optical depths > 3. Such high optical depths are rare: Heiles Sz Troland (2003) 


surveyed 79 21-cm line absorption spectra and, apart from two positions in the Galactic plane, only 


one exhibited an optical depth exceeding 3. Also, Stanimirovic et al. (2014) investigated whether 


the optical depth effect was significant in their inferring total H 2 content by comparing dust IR 
emission and 21-cm line emission in the Perseus molecular cloud by measuring 26 21-cm absorption 
line spectra in the immediate vicinity, and concluded that the optical depth effect was essentially 
negligible. 


4.4. Interpreting as dust property variations 


The apparently-low gas-to-dust ratios observed toward the isolated clouds in this study could 
be due to the amount of dust being overestimated, which could happen if the dust properties 
are different in the clouds than in the diffuse interstellar medium. Recall the gas-to-dust ratios 


are calculated assuming a value of the dust opacity, n, which we take from Li Sz Draine (2001). 


To explain a gas-to-dust ratio lower than the diffuse ISM by a factor of 3, we require that k 
is higher by a factor of 3. Variations of this amount were indeed inferred from the first sub 
millimeter measurements of a diffuse interstellar cloud and comparison to extinction from star 
counts (Bernard et al. 1999 Stepnik et al. 2003). A recent comparison of Herschel emission to 


near-infrared extinction for cold cores shows the ratio of 250 /rrn to 1.2 fj ,m optical depth is higher 


than in the diffuse ISM by a factor of 3 (Juvela et al. 2015); this is a ‘dust to dust’ comparison 
unaffected by ‘dark’ gas. The earlier subnnn and Planck results both indicate higher k, in the 
same direction and approximate amount required to explain the gas-to-dust variations derived in 
the present study. Observations of gas and dust in circumstellar disks yields k about a factor of 10 
lower (Beckwith et al.][l990); such a low k, if applied to the ISM, would yield gas-to-dust ratios 10 
times higher, which goes in the opposite direction from what is observed in the cloud cores and is 
inconsistent with the observed depletions of metals in the diffuse ISM. 

Because k is the ratio of the extinction cross section to the dust mass, it is sensitive to the 
size distribution of the large grains. For a compact (solid) particle, n is inversely proportional 
to the grain size, so grain growth would decrease k. But for coagulation of particles into fluffy 
aggregates, n apparently increases by a factor of 2-3 for silicate grains in the size range relevant to 


the ISM (Kohler et al. 2012). Variations in the dust properties are important to understand, not 


only because they may affect the inference of ‘dark gas’ presence, but also because of the role played 
by dust grains in the physics and chemistry of the interstellar medium, as the primary agents for 
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gas heating (through the photoelectric effect) and molecule formation (which occurs primarily on 
grain surfaces). 

Changes in dust composition are constrained by the availability of elements. Many heavy 
elements, including for example Fe and Si, are almost completely locked in grains (Dwek 1998), 
so unless the regions with different dust properties are enriched in heavy elements, any new grain 
construction would need to be done without the refractory material already depleted onto the grains. 
We have found no evidence of central sources in the clouds that could be responsible for heavy 
element enrichment. So for now we assume that regular cosmic abundances apply, and therefore 
that the changes in properties may be related to the grain’s geometry, degree of amorphization, or 
surface properties (like mantles). The carbonaceous and silicate grains may combine to form larger 
grains, for example, which may increase n. The wavelength-dependent emissivity may change as a 
result of any surface compositional change to the grains. This could occur when small grains cover 
the surfaces of the large grains. In this case it is not the ‘growth’ in size that matters but rather 
the change in chemical bonds affecting the manner in which energy is absorbed by and emitted 
from the surface. Carbonaceous material may evolve due to prolonged exposure to the interstellar 


radiation field, leading to a darkening over time (Cecchi-Pestellini et al. 2010, Jones et al. 2013). 


The nature of the silicate material may also be evolving in dark clouds. In order to explain the large 
values of /dark by increased grain emission, one would have to invoke the quite unlikely scenario 
of a greatly enhanced (factor ~ 3) metal abundance, all of which resides in the dust; this seems 
unlikely. 
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5. Conclusions 


We find strong variations in the ratio of gas to dust mass both within interstellar clouds and 
from cloud to cloud. Our work complements all-sky studies, which use large amounts of data but, 
when averaging together, give a misleading impression of uniformity, missing the true diversity of 
interstellar clouds. The three hypotheses that we present as potential explanations for the variations 
all have implications for the interstellar medium, and they also have observational tests that can be 
performed. The first two hypotheses are based on the idea that the gas column density may have 
been underestimated, while the true gas-to-dust mass ratio is constant everywhere. There would 


need to be extra gas, commonly referred to now as ‘dark gas’ following Grenier et al. (2005) not 
traced linearly by the 21-cm line emission. 

The possibility of molecular gas (first hypothesis) would not be a surprise given its prevalence 


in the diffuse interstellar medium as measured by UV absorption on random lines of sight (Rachford 


et al. 2002) and the widespread usage of dust images to trace the total gas column density even when 
in excess of the atomic gas column density (Planck Collaboration et al. 2011). For the brightest 


clouds in our sample, CO is detected and the H 2 /CO ratio appears ‘normal’ compared to galactic 
plane clouds. But for many clouds, and for all cloud peripheries, there is no CO and the upper 
limits show that any H 2 present is far in excess of that expected based on the CO emission. The 
implication would be that the clouds in this sample are in fact molecular clouds, despite being 
translucent and tenuous. The clouds are not gravitationally bound, with gravitational forces being 
exceeded by the combined influence of both thermal and non-thermal (turbulent) motions by factors 
of 10 2 to 10 3 . Observational tests can confirm or deny the presence of molecules, such as OH (which 
is readily formed wherever H 2 is present) in the clouds at abundance levels sufficient to explain the 
apparent ‘dark gas’ as H 2 . 


The second hypothesis is that the 21-cm line is optically thick, which is expected if the atomic 
gas is cold and has been proposed as an explanation for the apparent ‘dark gas’. We find that high 
optical depths are required, so that the 21-cm line is optically thick over entire interstellar clouds. 
The implications of this hypothesis are that 21-cm line observations are poor tracers of interstellar 
column density. The gas is likely to be cold and have optical depths > 3 in the 21-cm line over 
significant areas. This hypothesis can be tested with absorption line measurements toward radio 
sources. 


The third hypothesis is that the dust grain properties vary, in which case there may be no 
‘dark gas’ at all, but instead a new physical process that changes grain properties. If the grain 
properties evolve in translucent clouds, then the building blocks of molecular clouds also begin with 
these modified grain properties. Given the constraints of cosmic abundances and the high efficiency 
with which dust absorbs starlight already, changes in dust grain properties are unlikely to explain 
the observed gas-to-dust variations, but this hypothesis deserves to be explored and compared to 
independent observational data. 










34 - 


The Galactic Arecibo L-Band Feed Array HI (GALFA-HI) Survey data set was obtained with 
the Arecibo L-band Feed Array (ALFA) on the Arecibo 305m telescope. The Arecibo Observatory 
is operated by SRI International under a cooperative agreement with the National Science Foun¬ 
dation (AST-1100968), and in alliance with Ana G. Mndez-Universidad Metropolitana, and the 
Universities Space Research Association. 

Facilities: Planck, Arecibo, IRAS 
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